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Abstract 

We study the transition between sharp and smooth density distributions at 

the edges of Quantum Hall Liquids in the presence of interactions. We find 

that, for strong confining potentials, the edge of a v = 1 liquid is described 

by the Zp = 1 Fermi Liquid theory, even in the presence of interactions, a 

consequence of the chiral nature of the system. When the edge confining 

potential is decreased beyond a point, the edge undergoes a reconstruction 

and electrons start to deposit a distance ~ 2 magnetic lengths away from 

the initial QH Liquid. Within the Hartree-Fock approximation, a new pair 

of branches of gapless edge excitations is generated after the transition. We 

show that the transition is controlled by the balance between a long-ranged 

repulsive Hartree term and a short-ranged attractive exchange term. Such 

transition also occurs for Quantum Dots in the Quantum Hall Regime, and 

should be observable in resonant tunneling experiments. Electron tunneling 

into the reconstructed edge is also discussed. 
PACS:73.20.Dx 
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I. INTRODUCTION 



In Fractional Quantum Hall (FQH) states there are no bulk gapless excitations; the only 
gapless modes are edge states, which are responsible for non-trivial transport properties 
at low temperatures. Edge states arise naturally in real samples, as the two dimensional 
electron gas (2DEG) is confined in a finite region. The manner in which the 2DEG is confined 
determines the structure of the electronic density on the borders of the sample, and rich 
structures may appear. For smooth edges, strips of compressible and incompressible FQH 
states may be formed [fl]||. The structure of these smooth edges have been studied mainly by 
focusing on the electronic density distribution at large length scales, where it is reasonable 
to use a semi-classical approach . This approach consists in assuming that the electron 
density n varies slowly enough, so that one can use bulk values for the internal energy u(n), 
which has cusps for n corresponding to fractional filling factors. An improved calculation 
using Hartree-Fock approximation was done in Ref. || at finite temperature, which agrees 
very well with the electrostatic calculation |4]]. 

In the opposite extreme, the electron density varies sharply at the edge, and there is 
no room for the formation of incompressible strips. The sharp edge of the v = 1 state is 
described by a ID (chiral) Fermi Liquid, in which the occupation in the momentum space 
has a sharp drop at the Fermi momentum. This sharp drop in the momentum space is 
related to the fast drop of the electron density near the edge. 

One naturally questions how the sharp edge picture evolves into the smooth edge picture 
as the edge potential becomes smoother. As an interacting ID system, the Fermi edge of 
the v = 1 state may have the following possible singularities displayed in Figure 1, such as 
a Fermi Liquid singularity, with a Zp = 1 discontinuity (Fig. la) or a renormalized Zp < 1 
discontinuity (Fig. lb), or a Luttinger Liquid singularity (Fig. lc). One scenario is that the 
sharp edge and the smooth edge are connected by the distributions in Fig. lb or Fig. lc. As 
the edge potential becomes smoother, the occupation distributions in Fig. lb and Fig. lc 
also get smoother. According to this picture, the smooth edge of the v = 1 state contains one 
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branch of gapless edge excitations which is described by a renormalized Fermi Liquid. The 
smoother the edge, the stronger the renormalization. However, the calculations presented 
in this paper suggest a new scenario for some natural confining potentials. We find that the 
chiral nature of this one dimensional system play an important role in determining the form 
of the singularity. Due to the chirality, the Z F = 1 edge is very stable. Fig. la correctly 
describes the edge structure for a range of edge potentials even for interacting electrons. 
However, as the edge potential is smoothened beyond a certain point, the edge undergoes 
a reconstruction. A pair of edge branches moving in opposite directions is generated, and 
the occupation distribution in Fig. la changes into the one in Fig. 2, which contains three 
Fermi points. The occupation (n k ) has algebraic singularities at these Fermi points. This 
new scenario has received some support from exact calculations on small systems. 

In this paper we focus on the structure of the electronic occupation density distribution 
at the boundary of a v — 1 liquid in the presence of interactions. We work with a Hilbert 
space restricted to the first Landau level, and in the main part of the paper we assume that 
the spins are fully polarized. The droplet is confined by its interaction with an underlying 
positive background (one way to introduce a confining potential). The paper is organized as 
follows. In section II we introduce the ID interacting version for the problem, and discuss 
the importance of chirality in determining the structure of the edge singularity. In section 
III we present exact numerical results for small systems, which support the picture that, 
before a discontinuous transition occurs, the chiral edge system is reasonably well described 
within the Hartree-Fock approximation. In section IV we study the effects of this transition 
for Quantum Dots, including possible experimentally observable effects. Finally, in section 
V we investigate the consequences of the transition on the dynamics of edge excitations, and 
in section VI the consequences to electron tunneling into the reconstructed edges. 
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II. THE ID INTERACTING MODEL AND CONSEQUENCES OF CHIRALITY 



A system of interacting particles in a 2-D QH droplet can be mapped into a one dimen- 
sional problem by enumerating the single particle wavefunctions of the first Landau level. 
The Hamiltonian of the interacting theory is 



The dispersion ex,x' is determined by a background charge p(x), which we use to control the 
confining potential. The 0a's are the single particle wavefunctions, labeled by the quantum 
number A. For example, in the symmetric gauge, A stands for the angular momentum 
quantum number m, with <f) m (x,y) = ^-^le - ' 2 ^/ 2 and z = (throughout the paper 
we work in units of magnetic length = 1). The wavepacket <p m is centered in a circle of 
radius R = \/2m. In the Landau gauge, A denotes the linear momentum in the x direction 
k x , with cj) kx {x,y) = ^^i/ 2 e' fc3: e~^~ fe ^ ; l//2 , and the wavepacket <p kx is centered at y = k x 
(L is the size of a system subject to periodic boundary conditions). 

Let us consider for now backgrounds p(x) that are invariant under certain symmetry 
transformations, such as rotations (if we are studying a circular droplet, using the symmetric 
gauge) or translations along the x direction (if we are studying a long strip, using the Landau 
gauge). In this case we have ex,y = e A 5a,a'- Impurities break such symmetries, and their 
effect will be considered later in the paper. Because the interaction V(\xi — x 2 \) depends 
only on the distance between x\ and xi, it is also invariant under these symmetries, and 
thus we can rewrite the Hamiltonian as 



H = 2 £ A,A' c[c x , + ]T V XuX2! X 3M C Al C A2 C A3 C A4 , 



(1) 




where 




H = J2exc{c x + £ V(S, A, A') c 



: A+<5 c A c A' c A'+<5 



(2) 



A 5,A,A' 
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It is this ID interacting model that will be the basis of our study of the v = 1 droplet. The 
question we want to address is how to determine the ground state occupation number (c\c\) 
for this theory. 

Without loss of generality, let us focus now on the problem of a strip with length L 
and periodic boundary conditions (equivalently, a cylinder of circumference L), using the 
Landau gauge. The QH Fluid lies on the surface of the cylinder, between its left (L) and 
right (R) boundaries (see Fig. 3). 

In a typical ID interacting theory we have non-Fermi Liquid behavior; the Fermi discon- 
tinuity is destroyed by the interactions, and the system is better described as a Luttinger 
Liquid. Notice, however, that the Hamiltonian in Eq. (0) has a peculiar difference from 
the usual ID Hamiltonian of an interacting system in the sense that the effective scattering 
potential V depends not only in the momentum transfered 5, but also in the momentum 
configuration (i.e., A and A') of the scattered particles. The Luttinger Liquid behavior is 
caused by the coupling between particle-hole excitations in the two distinct Fermi points. 
Now, for our system described in Eq. (0), the two Fermi points correspond to \l and \r, 
at the two boundaries of our droplet. These two points are spatially separated, and the 
matrix elements for coupled particle-hole excitations near these points should go to zero as 
the distance between the boundaries is increased. In the limit of infinite separation, the two 
edges are decoupled, and we can describe the system as containing two different types of 
fermions, R and L, with one Fermi point each. More precisely, we can describe the particles 
as being in a Dirac sea that is filled as we move inwards to the bulk. Indeed, in topologies 
such as a simply connected droplet, like a disk, we only have one boundary, and the Dirac 
sea description is exact. Such theories fall within what we call "Chiral Luttinger Liquids". 

The ID chiral theory has a special property that the occupation distribution of the ground 
state can have a Fermi discontinuity (Fermi Liquid) even in the presence of interactions. In 
fact, for certain values of interaction strength and single particle dispersion, the ground 
state may have a perfect Fermi distribution with Zp = 1! This is because the momentum 
occupation with Zp = 1 is always an eigenstate of the interacting Hamiltonian (§) in the 
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limit of infinitely separated edges. This is easier seen in the filled Dirac sea description of 
the L and R fermions. Take, for example, the R branch, for which the unique minimum 
total momentum eigenstate is the one that has all single particle levels to the left of the edge 
occupied. Because total momentum commutes with the Hamiltonian, and this state is the 
only one with minimum total momentum, it must also be an eigenstate of energy, possibly 
the ground state for some edge potential. Notice that the occupation distribution of this 
state corresponds exactly to a Zp = 1 Fermi Liquid occupation. One should contrast this 
case with a non-chiral ID system, where clearly the Fermi gas distribution is an eigenstate 
of zero total momentum, but it is not the only one, and thus not necessarily an eigenstate 
of energy. Again, chirality plays a key role. 

The next step is to understand how the occupation distribution evolves as we smoothen 
the confining potential. One way to assemble a sharp distribution is by simply laying the 
electron gas on top of a similarly sharp positively charged background, and one way to try 
to destroy this sharp distribution is to smoothen the positive background. Notice that the 
perturbation we include by changing the background is not in the form of an additional 
interaction between the particles, but of a change in the one particle dispersion e\. We 
will show that the occupation distribution has the tendency to remain sharp, due to a 
balance between a repulsive long range Hartree term and an attractive short range Fock 
term, and also due to the special stability of the Zp = 1 chiral Fermi Liquid. The sharp 
distribution eventually becomes unstable, and the Fermi surface is destroyed, as a lump of 
particles detach and form two more edges which destroy the chirality. We will show that 
the Hartree-Fock approximation seems to contain the relevant ingredients to describe this 
transition. 

III. EXACT RESULTS FOR SMALL SYSTEMS 

In the Landau gauge, the dispersion due to the background charge, and the matrix 
elements V in Eq. (|2p for the Coulomb interaction are given, respectively, by: 
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e 2 roo poo e ~{y-k+y') 2 
£k = — / dy' p(y') / dy = In y 2 (3) 

6 J — oo J—oo y 7T 

and 

V(q, Ak) — — j- ~j= / dy e -^- Afe ) 2 / 2 2f (gy) (4) 

£ -t-' y Zll J — oo 

where e is the electron charge, e is the dielectric constant, q is the momentum transfer (5 
in Eq. (0)), Ak = k — k', p is the positive background density, and K is a modified Bessel 
function. Notice that e 2 /els (or e 2 /e, as we use units of Ib = 1) is the natural energy scale 
in the problem. In particular, the Hartree-Fock effective two body potential between two 
particles with momenta k\ and is given by 

V HF (h, k 2 ) = Vndh - k 2 \) + V^Qkt - k 2 \) (5) 

where the Hartree and exchange terms are obtained from Eq. (H) by setting q — > 0, AA: = 
k\ — k 2 and q = k\ — k 2 , Ak = respectively (with a factor of —1 for the exchange): 

1 P 2 f ,-(y-k 1 +k 2 ) 2 /2 

V H (\k 1 -k 2 \) = - — e - dy 6 - = lny 2 (6) 

LLi 6 J—oo y Z7T 

V ex (\h - k 2 \) = -—- e ( ^ ) K U — - — ) 

We have subtract a logarithmic divergence from the Hartree term (~ \nq\ q ^o), that is 
independent of ki and k 2 , and thus simply contributes to a constant in the energy. The 
Hartree contribution to the effective two body potential, Vh, is repulsive and long-ranged, 
whereas the one from exchange, V ex , is attractive and short-ranged. We will show that it is a 
balance between these two effective interactions that controls the short length scale behavior 
of the density distribution. 

We will proceed by first presenting exact numerical results for a small system, and then 
using these results to justify a picture that the short length scale behavior of the density 
distribution is controlled by the Hartree-Fock terms. 

We study the edge structure of a system that we divide into "edge" and "bulk" electrons 
(See Fig. 4). We consider just one edge, say, the R edge, and assume the "bulk" extends to 

7 




infinity in the opposite direction. The occupation of the "bulk" levels is fixed to be 1. Doing 
so, we can concentrate all the computations on the "edge" , as the effect of "bulk" electrons 
is simply reduced to a contribution to the one particle dispersion of the "edge" electrons. 
Such division presents no harm, as long as the edge excitations under consideration do not 
change the "bulk" occupation. 

Consider a strip geometry with L = 20 Ib and periodic boundary conditions (cylinder). 
The "edge" is composed of 10 electrons in 20 single particle states. We start with a sharp 
background, and smoothen it by changing the width w in which the density drops from the 
bulk value (p = l/2n, for v — 1) to zero (see Fig. 5). The "bulk" electrons contribute to 
an additional term in the dispersion. The effect of adding a variation in the positive charge 
density over a length scale w can be thought of as simply superimposing a dipole to the 
effective edge potential for a sharp edge, as shown in Fig. 5. 

Fig. 6 displays the energy levels for different total momentum K of the edge electrons 
(the sites, numbered from 1 to 20, are assigned k = to 19). For w < 8 l B the ground state 
has K = 45, i.e., all the electrons are packed up to one side and the edge is described by 
the Z F = 1 Fermi liquid. For w = 9 Ib the ground state is no longer the sharp configuration 
with minimum K = 45, but has moved to a configuration with K = 60 (see Fig. 7). The 
occupation number distribution is shown in Fig. 8 for w = 9 and 10 Ib- Notice the formation 
of a lump of electrons distant ~ 2 l B from the bulk (for L = 20, Ak = 1 corresponds to a 
distance 2n/L ~ 0.314 l B ). 

The exact calculation for a small system seems to indicate that the sharp edge is robust 
against the smoothening of the background charge, up to a point where there is a transition, 
and the density redistributes. We will show that the robustness of the sharp edge is a 
consequence of the attractive exchange, that tents to keep the edge packed. Eventually, as 
the strength of the dipole resulting from the smoothening of the background is increased, 
the short range attraction due to exchange can no longer sustain the edge sharp, and a lump 
of the electrons splits and forms a "puddle" near the minimum of the effective potential seen 
by the QH Liquid. Let's illustrate the point above by calculating the effective single particle 
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energy within the Hartree-Fock approximation for the distribution that has occupied levels 
for all negative momenta, i.e., a sharp R edge (we have centered the coordinate system on 
the edge). 



One should notice that the electronic occupation number n k = (c k c k ) differs from the elec- 
tronic charge density n(y) = (c* (y)c(y)) . The later is obtained from the former by a con- 
volution with a Gaussian of variance a 2 = 1/2 (the single level charge distribution). With 
this in mind, one can show that the background density that cancels the electronic charge 
density is exactly the one that makes e k + £#(&) = 0, as it should be expect. A sharp 
background charge distribution as shown in Fig. 5 makes a sharp electronic occupation 
distribution even more stable, as we have an extra dipole term, resulting from the difference 
between a sharp positive charge background and the electronic charge density of a sharp 
occupation distribution, and which favors the levels with negative k to remain occupied. 

As we change the background configuration, we alter the one particle dispersion e k . 
For the stability of this edge it is necessary that the effective single particle energy of any 
unoccupied level be higher than the one of any occupied level. Fig. 9 shows the effective 
single particle energy for unoccupied levels for different values of the parameter w, which 
measures the width which it takes the background to decrease from its bulk value to zero. 
The potential obtained for w = Ib is primarily due to the exchange term, which stabilizes 
a sharp edge. The exchange potential is short ranged, reaching zero within ~ 1.5 — 2 Ib] 
the overshoot for w = Ib is due to the dipole which results from the difference between 



e(k) = e k + £ H (fc) + ^ex(k) , 



(7) 



where 
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a sharp positive charge background and the electronic charge density of a sharp occupation 
distribution, as mentioned previously, whose contribution decays to zero as l/\k\ for large \k\. 
For w ~ 11 Ib the condition for stability is violated (the higher value for the w that marks 
the transition, as compared to the small system result, can be regarded as due to a finite size 
effect, to the Hartree-Fock approximation, or to both). It is then more advantageous to move 
electrons to the minimum e(k) locations. A simple picture is that particles start to escape 
from the sharp edge and start to deposit at a distance of order ~ 2 Ib away from the initial 
boundary. This separated lump brings in two new boundaries into the problem. These new 
boundaries break our previous chiral geometry, as we now have three Fermi points finitely 
separated. The three Fermi points describe two right-moving branches and one left-moving 
branch of edge excitations. From this point on one should expect that the interactions will 
destroy the Fermi Liquid singularities, and we will have three Luttinger singularities. 

Notice that we said the necessary condition for the stability of the sharp edge is that 
e(k) be larger for unoccupied states than for occupied states. But we have not yet argued it 
is sufficient. It is possible that even if this condition is satisfied one can have a ground state 
for the interacting problem different than the sharp edge, as the energy could be lowered 
by rearranging many particles. Worse, it is possible that the hopping terms, which couple 
different states in configuration space (and are not included within Hartree-Fock), would 
completely modify the picture. In particular, the Hartree-Fock approximation does not 
allow the distributions in Fig. lb and Fig. lc. Therefore, we cannot use the Hartree-Fock 
calculation alone to judge which of the distributions in Fig. 1 and Fig. 2 is realized after 
the transition. However, the exact diagonalization results for small systems that we have 
presented support the picture described in Fig. 2. This suggests the transition is mainly 
controlled by Hartree-Fock terms, and that the conclusions depicted from the single particle 
potential obtained within Hartree-Fock seem to be qualitatively correct. Certainly we cannot 
rule out the possibility that Fig. lb and Fig. lc might be realized for some other interaction 
and edge potential. 

To finalize this section, we present in Fig. 10 the spectrum and occupation numbers 
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calculated for the small system with w — 10 Is, but now within the Hartree-Fock approx- 
imation (the hopping or off-diagonal elements were suppressed). Compare the spectrum 
to the exact diagonalization for w = 101b displayed in Fig. 7. The occupation numbers 
in the Hartree-Fock approximation suggest the separation of part of the density from the 
main fluid, and the appearance of two more singularities. The hopping elements would take 
charge in redistributing the density, modifying the form of the singularity. 

IV. THE EDGE RECONSTRUCTION FOR QUANTUM DOTS 

The effect we describe in this paper is not particular to large systems. In fact, the 
exact results for small systems, which we used to support the Hartree-Fock picture, directly 
indicate that the transition occurs for finite systems. The edge density redistribution is a 
consequence of the balance between the confining potential, the repulsive Hartree term, and 
the short-ranged exchange, all of these present regardless of the size of the system. 

Quantum Dots are innately interesting systems for observing this transition. To begin 
with, because the number of electrons is small, the redistribution will involve a substantial 
part of the total number of particles in the dot, which can then be considered not simply 
an edge effect, but, in a way, a bulk effect as well. Secondly, experiments on resonant 
tunneling into Quantum Dots in the FQH regime should be sensitive to a transition involving 
a redistribution of the particle density, both because the energies of adding one electron to 
the dot on both sides of the transition should differ (which can be measured by the position 
of the resonant peaks), and because a change in the size of the dot will change the coupling 
to the probe leads as well. Third, the transition can be driven by altering the confining 
potential, either changing the voltage on a back-gate, or changing the magnetic field (which 
varies the radii of the orbits, and consequently the potential seen by each orbit). 

In this section we study some of the consequences of the transition as applied to Quantum 
Dots. We study systems with total number of particles up to N p = 70. We will use 
only the Hartree-Fock matrix elements, because both the system is not small enough for 
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exact diagonalization, and the Hartree-Fock approximation seemed to contain the essential 
elements to describe the transition. 

The energy eigenstates of the Hamiltonian in Eq. (Q) are also eigenstates of total mo- 
mentum. The true ground state is a superposition of different occupation number states in 
Fock space, all of them with the same total momentum. We will call hopping elements the 
terms in the Hamiltonian that couple different states in Fock space. The Hartree and Fock 
terms couple a configuration in Fock space to itself. Within Hartree-Fock, i.e., neglecting 
the hopping terms in the Hamiltonian, any occupation number state is an energy eigenstate. 
Finding the ground state is then equivalent to determining the configuration of particles 
that minimizes a classical energy function. Notice that within Hartree-Fock all (c^c^) are 
equal to either or 1. 

We focus on a disk geometry, which is more appropriate for describing a dot. The 
single particle states are labelled by the angular momentum quantum number. The matrix 
elements are given by 

Knl,m2,m3,m4 = (ml m3\V\m2 mi) , (8) 

where the state \m m') stands for a particle in the level labelled by m, and another in the 
level m! ( (z x ,z 2 \m m') = (j) m , m >(z x ,z 2 ) = j= jgaL e-^ ^ , with z Xy2 = Xl - 2 ^- 2 ). The 
matrix elements in the classical energy function that couple states m and m! are obtained 
from the Hartree and Fock terms: 

^mm 1 ^m,m,m' ,m' ^m' ,m,m,m' ■ (9) 

To obtain these coefficients it is easier to work in a basis in which V is diagonal. The 

z"z_ 



1 Z Z ' I 1 'i I 

\l,n) basis, in which (21,22^,71) = -7= t ~ e 2 ; where z± = , is such that 



(I, n\V\l', n') = V(l) 8i t u 5 n , n i. For the Coulomb interaction we have 

1 e 2 f dz\dz 2 _ 1 \z+\ 2n \z_\ 21 _ ]z+] 2 _| Z _| 



w 2el B J 7T 2 2\z-\ nll\ 
e 2 1 T(/ + 1/2) 

~ 7T B 4 r(z + i) ' 
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The Hartree and exchange terms are obtained, respectively, using 

(m m'\V\m m!) = ^ V(l) \{mm'\l,n)\ 2 (11) 

l,n 

and 

(m m\V\m m) = ^ (-1)' V(l) \(m m'\l,n)\ 2 . (12) 

l,n 

The confining potential is assumed to be parabolic, e(r) = |Ax 2 , with r the distance from 
the center of the dot, and k the strength of the confining potential. It is convenient to write 
k = ao-^lg 2 , so that «o is a dimensionless parameter, maintaining e 2 /eIb and as our, 
respectively, energy and length units. The single particle dispersion is given by 

1 e 2 rdz 2 „, l0 |z |2m 



e m — - o: 



f ^- 2\z\ 2 ^- e~\*\ 2 (13) 
J n ml 



2 el B 

e 2 r(m + 2) e 2 
e Lb r(m +1) e 

The energy function that must be minimized is 

£ = J]] Cm n m + J]] Kti,to' ^ m 5 (14) 

m m,m' 

where the or 1, constrained to J2n m = N p , the total number of particles. We 

obtained numerically V^^, for the first 80 levels (0 < m,m' < 79), and searched for the 
minimum of £ for different values of N p and a . We find that, depending on these parameters, 
the minimum energy configuration switches from a compacted to a separated droplet. 

In Fig. 11 we display the occupation of the orbits as function of «o for N p = 60 (occupied 
orbits are displayed in black, and unoccupied ones in white). For strong confining poten- 
tials the occupied levels are the ones with minimum angular momentum. As the confining 
strength is decreased, there is a transition, and unoccupied levels inside the dot appear. 
After the transition, hopping elements become important, and take charge in redistributing 
the occupation, which can then have partially filled levels. Fig. 12 displays the orbital occu- 
pation for fixed a , with N p varying from 70 to 30 electrons, where there is also a transition, 
with a separated droplet for smaller systems. 
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We would like to point out that the exchange term is of key importance in order to have 
a compacted dot solution. Notice that we have assumed that occupation is non-zero only 
for the first Landau level, and that the spins are fully polarized. One could argue that these 
assumptions alone can lead to a compacted drop, as a strong enough confining potential 
can always be chosen such that the ground state is the minimum total angular momentum 
solution even if one take only the repulsive Hartree term. This would be possible because 
the particles would be squeezed to the center, without being able to occupy higher Landau 
levels, or flip spin (see Fig. 13, where we repeat the calculation for Fig. 11 without the 
exchange term). In reality, as we increase the confining potential, there are two mechanisms 
which tend to lower the total energy of the dot, one by compacting the particles to the low 
angular momentum orbits in the first Landau level, with polarized spins, and the other by 
moving particles to a higher Landau level or opposite spin polarization state. By increasing 
the confining potential we enhance both of these effects. Therefore, in order to have a v — 1 
compacted dot we must have a window of a that allows totally packed dots with no higher 
Landau level occupation. It is here that exchange comes in play, providing an attractive 
interaction that lowers the bound on a to have a compacted dot, which opens that window. 

The ideas above can be expressed quantitatively. The lower bound on a , i.e., the min- 
imum confining strength necessary to keep an A^-particle dot compacted, can be obtained 
as follows. Within the Hartree approximation, ct™ m is obtained from the condition that 
the net electric field on the edge of the dot due to the electrons just balances the field due 
to the confining potential. The radial field due to the electrons diverges logarithmicaly, 
E r oc -^r \n(R/\ c ), where R is the radius of the droplet, and A c is an ultraviolet cut-off 
length scale. The field due to the confining potential is E r = —-jtcioR/Ib, so we find that 

el B 

a mm ^ jV" 1 / 2 \n( x y^ 2 ). Indeed, this dependence of a on N p fits very well the numerical 
results obtained when the exchange term is omitted (Hartree approximation), where we find 
a™ n ~ 0.118iV~ 1 / 2 ln(^) (see Fig. 14a). We find that a similar function dependance on 
N p reasonably fits the results obtained within Hartree-Fock, with a™ m ~ O.OSSiV^ 1 / 2 ln(^j-) 
for our range of N p (Fig. 14b). Notice that the attractive exchange term has the tendency 
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to keep the dot compated, lowering the value of a™ m . 

The upper bound on can be obtained by estimating the energy decrease of moving 
one particle from the edge to the center. The electrostatic energy of a disk of radius R 
and density po = 271 -V * s ^ disk = J^^(-^/^) 3 ' w hich gives an estimate for the electrostatic 
energy of adding one electron to the edge of the dot of 5S e d ge = ^^ L \[^~p- The electro- 
static energy cost of adding an electron to the center of the disk is 5£ center = ^\f2^W~ p . 
The total decrease in energy of moving one particle from the edge to the center is 
AS = (c*oN p — ^^V^y^Vp), which has to be < hu c (or < Szeeman) if we want to 
have occupation solely in the first Landau level (or with polarized spins). So we find 
max ^ 2=2. y/2 N~ 1/2 + nf c , (or %fnaa ) jV" 1 . The last term contains the ratio between 
the cyclotron (or Zeeman) and Coulomb energies. The condition for a compacted dot is 
a^ in < a < a^ ax . For GaAs -gf*- ~ OAy/B, with B the magnetic field in Tesla. So 
for reasonable values of B and N p in a dot, the term N~ 1 / 2 in a™ ax is the dominant one, 
which constrains the maximum possible N p to the one that makes a™ m ~ a™ ax , which 
gives i\T™ aa; ~ 106. If we use the value for a™ n given by the Hartree term alone we find 
jyrnax ^ Although these values are rough estimates, they should make it clear that the 
attractive exchange term plays a major role in opening up a window in «o for which there 
is a compacted dot solution. 

We have also performed the Hartree- Fock calculation for spin 1/2 electrons with Zeeman 
energy Szeeman ~^ + (simply to break the degeneracy between the two spin polarized con- 
figurations). For large confining potentials, both spin up and down electrons form compact 
droplets of filling fraction v = 1. However the droplet of, say, spin down electrons is smaller, 
and the electrons near the edge form a ferromagnetic state as pointed out in Ref. 0. As we 
decrease a , the separation between the spin-up edge and the spin-down edge increases. At 
even smaller oiq the spin-down electrons non-longer form a compact droplet which maybe a 
sign of FQH states. As «o decreases below a s ~ 0.53 N' 1 ^ 2 + 0.49 AT 1 , all the electrons are 
spin polarized (the coefficient in N~ 1 / 2 is approximatly the same obtained from the electro- 
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static consideration, and the one in iV" 1 shows the tendency of exchange to align spins as 
an effective Zeeman energy). Thus for a™ n < cto < a s , the electrons form a spin polarized 
compact droplet. The exchange term tends to align the spins, increasing further the size of 
the window of ao's such that the drop is compact. 

The next question is how to experimentally obtain a value for a that falls within the 
window above. In order to make the connection to real samples, we use the parabolic 
confining potential in Ref. @]. There they use V ex t{r) = ^m*u 2 r 2 , with m* the effective 
electron mass in GaAs, and Tiuj = 1.6meV for the particular device. The relation between 

2 

our dimensionless ao (or alternatively, given in units of = 1 and Z# = 1) to this confining 
potential is obtained by equating cuo^l^ 2 to m*uj 2 , which gives a = ctt^W^ k , or a ~ 
0.376 B~ 3 / 2 , B in Tesla. For N p = 40, for example, we find that B = 2.5T will yield a value 
of «o in the allowed window, so that the dot occupation will be the one of a compacted 
v = 1 droplet. As we increase B beyond 3. IT the edge will undergo a reconstruction. 

We now turn into the possibility of probing experimentally the transition between a 
compacted and an expanded dot. In order to make a clear connection between this expansion 
effect and experimentally observable quantities, we describe below the implications of the 
effect to tunneling experiments into Quantum Dots. In resonant tunneling experiments, the 
energy difference between the ground states of an iV + 1 and N electron system, [i(N), can 
be probed by tunneling in and out of the dot a single electron at a time, when the Fermi level 



of the electrodes become resonant with the quantum level of the dot [pHlO|j. By following 
a peak, the dependence of the chemical potential on the magnetic field can be observed. 
In Fig. 15 we calculated, within the Hartree-Fock approximation, this dependence of the 
chemical potential (in meV) on B (in T) for dots with N p from 35 to 38 . The "sawtooth" 
for B < 2.5T corresponds to spin down electrons being flipped and moved from the center 
to the edge of the dot. The region between roughly 2.5T and 3T is the window of magnetic 
fields for which the electrons form a compact v = 1 droplet. The "dislocation" near B = 3T 
corresponds to the reconstruction of the electron number occupation, marking the transition 
from the compacted to the expanded configuration. 
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In addition to this anomaly in the peak position vs. B, the expansion of the size of 
the dot will also increase its coupling to the probing leads, as this coupling depends on 
the distance between leads and island. The tunneling current should then increase for an 
expanded dot. In Fig. 16 we show the dependence of the size of the droplet (measured as 
the radius of the orbit of the outmost electron) on a and on the number of particles N p in 
the dot. In Fig. 16a we show the size of a 60 electron droplet as a function of ceo, and in 
Fig. 16b we fixed the value of a and varied the number of electrons from 40 to 70. The 
effect on the amplitute of the resonant peak, together with the anomaly in the peak position, 
should be a signature that a transition is indeed occurring in the occupation density of the 
Quantum Dot. 



V. THE EDGE MODES AFTER THE TRANSITION 

As we have seen in the previous sections, after smoothening the edge potential enough, a 
transition takes place, and the Fermi Liquid occupation density gives way to a more complex 
state. In terms of the electronic occupation distribution, the new state looks as if electrons 
start to deposit a certain distance away from the bulk of the QH Liquid. The QH "puddle" 
that is formed, as mentioned above, brings in two more boundaries for each edge, and we 
then have three singularities. The interactions take charge in destroying the Fermi Liquid 
discontinuity, as the three singularities are finitely separated. 

These three singularities can be related to three branches of gapless modes, which corre- 
spond to particle-hole excitations near each of the singularities. An intuitive way to visual- 
ize the three branches is by considering the occupation number after the transition simply 
within Hartree-Fock, which would look as in Fig. 17. There, the three edges are clearly 
identified. Correlations destroy the Fermi discontinuities, but we will still have the three 
Luttinger Liquid singularities at the three Fermi points. Notice that the pair of branches 
that is added always have opposite chirality. An edge that had one right-moving branch 
before the transition, for example, will have two right-moving branches and one left-moving 
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branch. All these three branches are strongly coupled. One clear experimental consequence 
of now having one branch moving in the opposite direction (the one left-moving branch 
in the originally right-moving edge, for example) is that one could probe such excitations, 
where originally there was none. We will show, however, that the presence of impurities 
localize two of the three branches, and such back-propagating modes cannot be observed 
beyond the localization length. 

We will use the bosonized description of the edge states in the FQH regime presented in 
Ref. [Ill, and which we summarize below for our particular case of v = 1. Let 4>r,l be two 



fields, described by the Lagrangian density 

£r,l = ~r 9 x 4>r,l (±d t - vd x )<p R)L (15) 

47T 

(v is the velocity of the excitations) and the equal-time commutation relations 

[(pR,L(t,x) , <j) RiL (t,y)] = ±wr sgn(x - y) . (16) 

Left and right moving electron operators can be written as ^r : l(x, t) =: e^ R ' h ^ x,t ' :, which 
can be shown to satisfy the correct anticommutation relations [JHJ]. The electron density is 
given by pr^l = <9x0r,l, and the Hamiltonian is 



H - V 



/ dx p\ L . (17) 



Consider three edge branches as depicted in Fig. 18, labelled by % = 1,2,3. We can 
generalize the description above to include several branches, writing the following Lagrangian 
density: 

£ = \ K ii d t<Pid x <pj - Vij d^d^j] , (18) 

where the matrices K and V contain, respectively, information on the direction of propaga- 
tion (chirality) of each branch, and interactions between the branches (including a diagonal 
term containing the velocities). For this analysis let us assume two R branches (1 and 3) 
and one L branch (2), such that 



K 



( 1 0\ 
0-10 

Vo o ly 



(19) 



Electron operators can be written as 



m L oc e^i'*** 



(20) 
(21) 

where the Lj's are integers satisfying = 1. The Lagrangian in (fl8|) is not the most 

general one. We left out four fermion terms that cannot be written as the product of two 
densities. We know that it is the interactions between the densities that are responsible for 
changing the Fermi edge singularity, and here we will concentrate on this particular effect of 
interactions. We have also assumed a local interaction in the densities. However, the long 
range interaction can be easily included by allowing Vij to have momentum dependence. In 
particular, the long range Coulomb interaction contributes to a term J2k ^k(J2i Pk,i)(J2i P-k,i), 
where A& oc In A; is the Fourier transformation of the 1/r Coulomb interaction and p^j the 
Fourier component of the density of the i th branch. We see at long distances the most 
important interaction term is the one that involves only the total charge density J2 Pi- This 
contribution is 

2 



\ i / 1,3 



(22) 



which, when summed to the velocity terms, gives the total V matrix 

/^i + A A A \ 

V = 



(23) 



A v 2 + A A 
V A A v 3 + \J 
Here, for simplicity, we have assumed that A is a large constant independent of momentum. 
This will be the case if the Coulomb interaction is screened at a long distance (e.g. by gates 
nearby). We will focus primarily in the case where A 3> f's, i.e., strongly coupled branches. 
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Also, if the system has particle-hole symmetry, then v\ = V3. We start by rewriting the 
Lagrangian (|TJ) in terms of new fields fa = J2j Uy <pj that simultaneously diagonalize K 
and V. Furthermore, we would like to keep, for convenience, 

K = {U T Y l KU- 1 = K (24) 

so that the commutation relations of the fas are the same as the ones for the fas. The 
transformation matrix U for A ^> v's is 





t 


1 


1 


1 \ 






1 

V~2 




1 




\ 


1 

V2 








The mode in the first line of U, 4>\ = + 4>2 + 03, is simply the total charge density mode 
(equivalently, pi = pi + P2 + P3, as p = d x fa). What we show next is that, once we add to 
the Lagrangian (|i~8D electron scattering terms due to the presence of impurities, this total 
charge mode is left unperturbed, and the other two will localize. 

The scattering terms between electron operators in the three edges that can be added 
to the Lagrangian are bosonic couplings with zero charge. These can be written as Tl = 
e % ^Zi l i4>i — e l Ylij L ^ K t]<t>i ^ w ]j ere now the L { 7 s are integers satisfying J2i^i = (Tl is bosonic 
and neutral). In terms of the rotated fields fas, Tl = e*^*^, where k = Y^jljUZ 1 . Let 
us focus on the most relevant TVs; the naive RG dimension can be obtained from the Tl 
correlations 

(Tl(t = 0, x) T L (t = 0, x = 0)> oc e~^ Ui ^(t=o,x) &(t=o,*=o)> (2g) 

OC X~^i li = X~ lL , 

Writing 7^ in terms of the L's we obtain: 

Tl = E? = ^~ (27) 

= i T u- 1 (u~ l ) T i 

= L T K\J- X (\J- 1 ) T KL 
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= L T U T K 2 UL 
= L T (U T U)L , 

where we used K 2 = K 2 = I. It is easy to show that the minimum 7^, with J2i Li = 0, are 
given by: 

/ 1 \ / \ 



L = ± 



-1 
\ / 



and L — ± 



-1 
x 1 J 



(28) 



which correspond to T L operators that transfer charges between the center branch (L) to 
the two side branches (R). In terms of the Z's, we have: 

/0\ /ON 



Z = (U- l ) T KL = KUL = ± 



1 



v-73/ 



and ± 



1 



The Hamiltonian density with these most relevant Tl terms added is 



(29) 



7i = 7~Co + Tix 



(30) 



where 



and 



H = £ (d^) 2 + ^ (^0 2 ) 2 + ^ (^0 3 ) 2 



47T 



47T 



47T 



• 4*2 ■ < f > 2~4 > 3 

H T = i+{x)e + £_(x)e V2 +h.c. 



(31) 



(32) 



The £± describe the random tunneling coupling due to impurities, with correlations 
(e±(^±(y)> = A±(5(x-y). 

Notice that 0i remains free, as the added tunneling terms do not depend on it. This one 
component (total charge, as fi\ = p\ + p 2 + Pz) behaves just like the one branch before the 
transition. The other two, we will argue below, should be localized because of the impurities. 
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Let 



hi ± ^3) / y/2, which obey the commutation relations 



*>±(t,x) , (f)±(t,y) 
l±(t,x) , <j> T (t,y) 



(33) 



-in sgn(x — y) 



We can identify fi± = ^ as the conjugate momenta to We can rewrite the Hamiltonian 
for 2 and 03 in terms of 0±: 



2,3 



47T 



(«9 X + ) 2 + (^0_) 2 + £+(z)e^+ + £-(z)e^- + ff.c. 



(34) 



where we assume that V2 ~ W3 ~ u. This is a more complicated version of a Sine-Gordon 
(SG) Hamiltonian density with position dependent coupling, as it involves self-interactions 
in both a field and its conjugate momentum. 

If we had only one of £ + or we would have a Hamiltonian for a simple SG with 
position dependent coupling, which we could write as 

v 



n 



4-7T 



(27rn) 2 + (d x <PY + £(x)e l * + H.c. 



Working in units of v — 1, and rescaling the fields as n' 
keep the commutation relations unchanged), we have 

1 



27rn and 



H =2 



n ,2 + (^0') 2 +i(x)e 1 ^' + H.c. 



(35) 
0/\/27r (which 

(36) 



with g = \[2n. This problem, equivalent to a Coulomb gas with position dependent chemical 
potential, was studied in refs. fl3] , |l4]| . The impurity coupling A± is relevant for g < \fQn 
(in the constant coupling constant or chemical potential Coulomb gas, the condition is 
g < \/87r). This is indeed our case, and therefore the presence of impurities localizes the 
other two branches of excitations represented by 02 and 03. 

Notice that what we have done above is equivalent to understanding the RG flows in 
the planes A + = and A_ = 0, and this implies A± are relevant in all directions around 
A± = if g < \fQn. The RG flows to a strong fixed point when both A± 7^ 0. It is possible 
that this strong fixed point is a localized state, motivated by the flow when A_ = 0. The 
properties of this strong fixed point will be the subject of further studies. 
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VI. TUNNELLING INTO RECONSTRUCTED EDGES 



The formalism in the last section can be used to study the electron tunneling into the 
reconstructed edges. The electron propagator in time in general has a form 

< ct(t)c(0) >oc -L 

with ji given in Eq. (P7|) . But now the Lj's satisfy J2i Li = 1. In the limit X/v ^> 1 U is 
given by Eq. The minimum exponent is Jl — 1 (the Fermi liquid value) for electrons 

described by L = (1,-1,1). This configuration corresponds to adding two electrons on 
the two side branches and removing one electron from the center branch. Adding a single 
electron to the side branch leads to a exponent 7l = 2 and to the center branch 7l = 3. 

Let us consider tunneling between two reconstructed edges. At very low temperatures 
and low voltages, the electron with the configuration L = (1, —1, 1) will dominate the 
tunneling and leads to a linear I — V curve [15f|, since / cx V 2lL ~ l and 7l = 1. At higher 



voltages, depending on the sample geometry, it may be easier for an electron to just tunnel 
into the side branch (with configuration L = (1, 0, 0)). In this case I oc V 3 and (dl /dV)v=o oc 
T 2 . 

The above discussion also apply the reconstructed edges of Laughlin states of filling 
fraction 1/m. But now for the L = (1,-1,1) electron the exponent jl = m. 7l = 2m 
for L = (1,0,0) and 7^ = 3m for L = (0,1,0). We see that in the limit X/v ^> 1 the 
minimum exponent in the electron propagator is not affected by the edge reconstruction. 
This result is valid even when more then one pair of edge branches are generated. This is 
because adding electrons of configuration L = (1, — 1, 1, — 1, 1) just displace all the edge 
branches by the same amount. Thus the electron of L = (1,— 1,1,— 1,...,1) just couples 
to the total density Pi an d do not couple to other neutral modes. In the limit X/v 3> 1, 
the total density mode decouples from other neutral modes. This is the reason why the 
L = (1, —1, 1, —1, 1) electron always have the exponent jl = th. We would like to stress 
that the above result is valid only at low energies (energies below the smallest Fermi energy 
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of generated edge branches). The high energy behavior of the electron propagator is not 
clear. In that case it is probably better to view the edge region as a compressible gas. 

For tunneling between two reconstructed FQH edges, we expect / oc \Z 2m_1 and 
(dl /dV)v=o oc T 2m at low voltages and low temperatures. This is consistent with a re- 
cent experiment on tunneling between (smooth) edges of 1/3 FQH states. []nj] 

VII. CONCLUSION 

In this paper we studied the electronic density of Quantum Hall Liquids, focusing on 
short length scales that are comparable with the magnetic length. We found that sharp 
electronic occupation densities, corresponding to a Zp = 1 Fermi Liquid, is possible because 
of the chiral nature of the system. This sharp distribution is stable against variations in 
the confining potential up to a certain point, beyond which it undergoes a transition and 
electrons start to separate from the bulk and deposit a distance ~ 2 lg away. The transition is 
shown to be qualitatively described within the Hartree-Fock approximation. The separation 
generates a pair of branches of edge states that move in opposite directions. 

We would like to remark that the separated electrons do not form any fractional quantum 
Hall state. This is because the separation between the Fermi edges is always of order 
magnetic length for realistic potentials. In this case 1/3, 1/5,... states are all described by 
Luttinger liquid and are indistinguishable. 

We also present results for Quantum Dots, where we predict that this effect of edge 
separation can be related to an expansion of the dot, which could be experimentally observed. 

The authors would like to thank Dmitrii Chklovskii, Yong Baek Kim, David Abusch, 
Olivier Klein, Paul Belk, Marc Kastner, and Ray Ashoori for useful discussions. This work 
is supported by the NSF Grant No. DMR-91-14553. 
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FIGURE CAPTIONS 



Figure 1. Possible singularities for the momentum occupation distribution of a v — 1 state: 
(a) Zp — 1 Fermi Liquid singularity, (b) Zp < 1 Fermi Liquid singularity, and Luttinger 
Liquid singularity. 

Figure 2. Momentum occupation distribution after reconstruction, with three singulari- 
ties. The addition of two singularities correspond to the addition of two branches of opposite 
moving edge excitations. 

Figure 3. Cylindrical geometry, equivalent to a strip of length L and periodic boundary 
conditions, where it is convenient to use the Landau Gauge. The QH Liquid (shaded area) 
lies on the surface of the cylinder, between its left and right edges at and \ R . 

Figure 4. The v — 1 droplet is divide into "bulk" , where all the states are fully occupied, 
and "edge" , in which states can be partially occupied. The division allows one to focus the 
computations solely on the "edge" electrons, with the "bulk" simply contributing to the one 
particle dispersion. This sort of division is not unique, as one can adjust the position of the 
boundary between the two regions; this boundary can be moved as long as the sites on the 
"edge" side of it are all fully occupied. 

Figure 5. Smoothened background charge density, which over a width w drops from its 
bulk value to zero. Such density can be written as the superposition of a sharp density 
profile tp a dipole term, which is used to tune the confining potential as function of w. 

Figure 6. Energy eigenstates obtained as function of the total momentum K of the edge 
electrons for values of w ranging from to 7 Ib- Notice that the ground state for this range 
has = K = 45, the minimum momentum configuration for 10 electrons (the 20 sites used in 
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the calculation are assigned k values from to 19). Also, notice that the energy levels of 
states of higher K are pulled down as w increases. 

Figure 7. For w > 8 Ib the ground state configuration is no longer the sharp edge 
distribution. The transition occurs near w = 8 Ib, and for w = 9 Ib is already fully 
developed, with the ground state momentum moved to K = 60. 

Figure 8. Occupation numbers for the ground state past the transition, for w = 9 and 
10 Ib- Notice how a lump of density moves away from the main body of QH fluid. This 
density profile (which goes from its bulk value to zero by first decreasing, then increasing, 
and finally decreasing again) is the signature of the existence of now three singularities, and 
thus three branches of edge excitations. 

Figure 9. The effective single particle potential calculated within the Hartree-Fock ap- 
proximation for different widths w, decreasing from Ib to 15 Ib in steps of 1 Ib- Notice 
that starting at w ~ 11 Ib, the condition for stability of a sharp edge is violated, as there are 
locations with smaller effective potential than the one at the edge of the sharp occupation 
density. 

Figure 10. Energy levels and occupation numbers for w = 101b calculated within the 
Hartree-Fock approximation (the hoping or off-diagonal elements were suppressed). 

Figure 11. Occupation number of the angular momentum states as a function of «o 
for N p = 60, calculated within the Hartree-Fock approximation. The occupied orbits are 
represented in black, and the unoccupied ones are shown in white. 

Figure 12. Occupation number of the angular momentum states as a function of N p for 
«o = 6.25 x 10~ 2 , calculated within the Hartree-Fock approximation. 
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Figure 13. Occupation number of the angular momentum states as a function of cto for 
N p = 60, calculated within the Hartree approximation. 

Figure 14. Minimum «o necessary to keep an A^-particle droplet compacted, calculated 
using the Hartree (a) and Hartree-Fock (b) approximation. The solid line is the best curve 
fit consistent with an electrostatic (Hartree) argument. 

Figure 15. Dependence on the magnetic field B of the energy cost to add one more 
particle, /i, to an island with 35 (lowest curve), 36, 37 and 38 (highest curve) electrons. 
The "sawtooth" corresponds to 1 < v < 2, where electrons are spin flipped and taken from 
the center to the edge as the magnetic field is increased. The "dislocations" near B ~ 3T 
correspond to the transition between compacted and separated dots. The region in between 
(B ~ 2.5 to 3T) is the window for which the dot is a compact v — 1 droplet. These results 
were obtained for the parabolic confining potential of the sample studied in Ref. [7J], with 
a = 0.3765- 3 / 2 (B in Tesla). 

Figure 16. Radius of the QH droplet for (a) fixed N p and varying a 0) an d (b) fixed «o 
and varying N p . One can vary ao by changing the strength of the confining potential, or by 
changing the magnetic field. For the parabolic confining potential of the sample studied in 
Ref. 0, the parameter a = 0.376B~ 3 / 2 (B in Tesla). 

Figure 17. Momentum occupation distribution after reconstruction, if calculated only 
within the Hartree-Fock approximation. The addition of two singularities correspond to the 
addition of two branches of opposite moving edge excitations. 

Figure 18. Three branches of edge excitations, two right-moving (1 and 3, on the sides) 
and one left-moving (2, in the center). 
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